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We report on observation of ferroelastic domain structure in single crystals of multiferroic 
CaMn70i2 at room temperature. Two types of ferroelastic domain wall are found, consistent with 
the material’s rhombohedral symmetry that is reduced from cubic symmetry at higher tempera¬ 
tures. Using Raman spectroscopy along with other measurements, we develop a systematic method 
to determine the microscopic domain orientation. Moreover, we hnd a switching behavior of the 
domains, which allows us to detwin the crystals conveniently at room temperature using a mod¬ 
erate uniaxial compression. Our result paves the way for further spectroscopic study and domain 
engineering in CaMn70i2. 

PACS numbers: 75.85.+t, 77.80.Dj, 78.30.-j 


I. INTRODUCTION 

Ferroic (ferroelectric, ferromagnetic, ferroelastic) do¬ 
main walls have aroused persistent research interest due 
to their significance in both fundamental research and 
promising applicationsIn so-called multiferroic ma¬ 
terials, different ferroic order parameters coexist and ex¬ 
hibit mutual coupling, hence allowing for the manipula¬ 
tion of one ferroic property through another. This co¬ 
existence gives rise to composite domain walls^“— that 
might be a key to utilizing the mutual controllability of 
ferroic properties in applications. Furthermore, domain 
walls can exhibit distinctly different properties from the 
bulk^“— leading to the possibility of using the domain 
walls as device. In addition to the intensively studied 
system BiFeOs^ intriguing domain structures have been 
observed and tuned in manganitesi^i^iii”— 

CaMn70i2 is a “type-II” multiferroic materiaU^ with 
very large ferroelectric polarization induced by mag¬ 
netic order, and has been the subject of considerable 
recent research effortUnlike in many widely stud¬ 
ied ferroelectric materials, the occurrence of ferroelec- 
tricity in CaMn70i2 is preceded by a ferroelastic struc¬ 
tural phase transition at higher temperatures, making 
it an ideal platform to study the ferroelastic and ferro¬ 
electric domain structures separately. At high tempera¬ 
tures, CaMn70i2 possesses the AC3B4O12 cubic struc¬ 
ture which is a derivative of simple perovskite ABOs^^ 
Upon cooling, a first-order phase transition occurs at 
Ts 440 K. The four Mn^-^^T ions in each formula unit 
are charge-ordered into three Mn^+ and one Mn^+ ions. 
The body diagonal of the high-temperature cubic cell 
that runs through the Mn^+ ions shrinks a little bit, and 
becomes the c-axis of the new rhombohedral unit cell in 
hexagonal basis This hexagonal c-axis plays an im¬ 
portant role that not only determines the direction of the 
incommensurate orbital order established below Tq = 250 
but also sets the direction of the giant improper fer¬ 
roelectric polarization that arises from the helical mag¬ 
netic order below Tati = 90 The pseudo-cubic 


cell below Tg, along with the shortened body diagonal, 
are denoted by yellow/blue cubes and magenta lines, re¬ 
spectively, in Figs.[Ha-b), omitting most details. 

Since current methods of growing single crystals of 
CaMn70i2 all take place at temperatures far above the 
ferroelastic structural transition temperature 
one can reasonably expect ferroelastic domains to form 
upon cooling the crystals to room temperature, as a 
result of simultaneous nucleation from different parts 
in the crystal, just like those in YBa2Cu306+(?^ and 
BaFe2As2^^ In this paper, we focus on domain struc¬ 
tures in CaMn70i2 at room temperature, where the com¬ 
pound is rhombohedral but not ferroelectric. We present 
direct observations of two types of stripe-like ferroelastic 
domain structures in single crystals. By using polarized- 
light microscopy, Raman spectroscopy and stylus surface 
profiler, we study the ferroelastic domain structures we 
observe. We propose a method to uniquely determine 
surface domain structures using Raman spectroscopy in 
combination with stylus profiler. Moreover, we identify a 
switching behavior of the ferroelastic domains, and show 
that the crystals can be mechanically detwinned by a 
moderate uniaxial compression. These findings are im¬ 
portant for further spectroscopic studies of CaMn70i2 
that require single-domain samples, and may facilitate 
future investigations of domain and domain-wall proper¬ 
ties in multifunctional oxides. 


II. EXPERIMENTAL METHODS 

High-quality cube-shaped single crystals of CaMn70i2 
were grown with a flux-reaction method^ at a cooling 
rate of 5 °C/h. Natural facets of the crystals are par¬ 
allel to crystallographic { 100 }^. planes, where the sub¬ 
script “c” denotes pseudo-cubic notation. The samples 
were characterized as described elsewhereJ^ Optical im¬ 
ages of crystal surfaces were taken with a polarized-light 
microscope Olympus BX 5 I with polarizer and analyzer 
set in an almost perpendicular configuration. Differential 
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interference contrast apparatus was installed to enhance 
image contrast. Surface profiles of crystal facets were 
measured with a KLA-Tencor P-6 Stylus Profiler using a 
contact force of 0.5 mg. Raman scattering measurements 
were performed in a back-scattering confocal geometry 
using the 632.8 nm line of a He-Ne laser for excitation. 
The diameter of the focused laser spot is estimated to 
be less than 5 microns. A Horiba Jobin Yvon Lab RAM 
HR Evolution spectrometer, equipped with a 600 gr/mm 
grating and a liquid-nitrogen-cooled CCD detector, was 
used to analyze the Raman spectra. 


III. RESULTS AND DISCUSSIONS 
A. Possible domain-wall orientations 

Figures [Hc-d) display typical polarized-light optical 
images of multi-domain samples. Ferroelastic domains 
manifest themselves in a regular bright-and-dark stripe 
pattern. While only the crystals’ top faces are shown, the 
domain structures actually extend to span the entire crys¬ 
tals, z.e., a consistent pattern is found on the side faces as 
well. This is illustrated (with exaggerated rhombohedral 
distortion of the domains) in Fig. [TJe) and Fig. [Hf), in 
which the slab-like domains are found to stack along the 
pseudo-cubic (lOO)c and (llO)c directions, respectively. 
On slightly twinned samples, especially those grown at a 
cooling rate slower than 5 °C/h, we can occasionally find 
stripes near the edges of the faces that do not span the 
entire crystal. 

When the crystals are heated to temperatures above 
Ts, the stripe patterns gradually disappear, consistent 
with the recovery of a single-domain cubic structure 
above Tg. Cycling the temperature through Tg can each 
time lead to a completely different stripe pattern in rhom¬ 
bohedral phase, which indicates that the domain forma¬ 
tion is not pinned by disorder or defects, and it in turn 
confirms the high quality of our samples. On the other 
hand, when the samples are cooled down to cryogenic 
temperatures, the stripe pattern remains even below T/vi 
= 90 K in the ferroelectric phase. Hence it is possible to 
control low-temperature ferroelectric domain structures 
by pre-setting a desired ferroelastic domain structure in 
the paraelectric phase. 

Grown at temperatures far above Tg and then cooled 
to room temperature, most crystals contain multiple fer¬ 
roelastic domains due to simultaneous rhombohedral dis¬ 
tortions that nucleates from different parts of the sam¬ 
ples. Even though as-grown crystals can occasionally 
be found in a single-domain state^d^ further manipula¬ 
tions may affect the domain structure. For example, to 
acquire pure Ag Raman spectra in CaMn 70 i 2 , the prepa¬ 
ration of a polished surface is required^ and the polish¬ 
ing process will inevitably exert mechanical stress onto 
the sample, thus raising the risk of twinning it. We will 
discuss domain-switching behavior under external forces 
and how to utilize it to detwin crystals later. 
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FIG. 1. (Color online) (a), (b) Schematics of {100}c and 
{110}c type domain walls in CaMn70i2. The blue and yellow 
cubes indicate different domains with the magenta lines indi¬ 
cating the shortened body diagonals. Domain walls are de¬ 
noted by the planes separating the cubes, (c), (d) Polarized- 
light optical images of single crystals with {100}c and {110}c 
domain walls. In addition to the domain structures which ap¬ 
pear as regular bright-and-dark stripe patterns, growth ter¬ 
races are seen especially in (c), but they do not seem to affect 
the domain distribution, (e), (f) Three-dimensional illustra¬ 
tions of the crystals in (c) and (d), respectively, with blue 
and yellow slabs indicating the two different domains. Sur¬ 
faces with four inequivalent domain structures are labeled as 
A, B, C, and D. 

Domain-wall orientations in ferroelastic materials can 
be understood by the equilibrium boundary condition 
(strain compatibility), which can be written as the 
following 
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Y =0’ (1) 

where S and S' are the strain tensors of two adjacent 
domains. Indices i and j = 1, 2, and 3 denote Cartesian 
coordinates, and all possible (xi,X 2 ,X 3 ) that satisfy the 
equation constitute permissible boundaries between the 
two domains. 

The rhombohedral phase of CaMn 70 i 2 belongs to 
the ferroelastic species m3F3 with four possible domain 
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variants^ which is the result of symmetry lowering from 
Im3 to R3i2i The four different spontaneous strain ten¬ 
sors can be written in the form>^ 
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with the corresponding shortened body diagonals lying 
in direction [lll]c, and respectively. 

When any two of these four domain variants meet, the 
solutions to Eq. [T] correspond to one of the equivalent 
crystallographic planes {100}c and {110}c^ which is in 
perfect agreement with our observations in Fig. [H The 
two different types of domain walls are illustrated in 
Figs.^Ja-b). Similar domain structures have been found 
in BiFeOs^ BaTiOs^ LaAlOs^ etc. All of these per- 
ovskites have rhombohedrally distorted phase, the strain 
tensors of which are in the same forn>2i as in our case. 
Additional domain variants may exist when ferroelectric- 
ity sets in, but the number of ferroelastic variants will 
remain to be four since 180° ferroelectric domain walls 
are not ferroelastic 

There are a total of four inequivalent surface domain 
structures on any {100}c face of a cube-shaped crystal. 
They are labeled as A, B, C, and D in the illustrations in 
Figs.[He-f). On an A-type face one finds (lOO)c stripes in 
the optical image, which continue onto adjacent faces as 
(100)c stripes; a B-type face exhibits no stripes since it is 
single-domain; a C-type face has (llO)c stripes; a D-type 
face has (lOO)c stripes, but unlike the A-type, the stripes 
continue onto adjacent faces as (llO)c stripes. Therefore, 
for a given crystal with sufficiently large single-domain 
volumes that span the entire crystal, one can obtain a 
rough idea about its domain structure by visually in¬ 
specting all of its six faces under a polarized-light micro¬ 
scope. 


B. Identification of individual domain orientation 

Optical images of crystal faces can only provide infor¬ 
mation about the relative orientations of the hexagonal 
c-axis in adjacent domains (Figs.[Ha-b)). Other methods 
are required to reveal the absolute orientation of individ¬ 
ual domains. We reported in our previous worki^ that 
Raman scattering can detect Ag and Eg optical phonons 
separately with parallel and perpendicular combinations 
of incoming- and scattered-photon polarizations, respec¬ 
tively, when at least one of the polarizations is parallel to 
the hexagonal c-axis. If none of the polarizations is along 
the hexagonal c-axis, the acquired Raman spectrum will 
be a weighted combination of Ag and Eg signals The 
relative intensities of Ag and Eg phonon peaks depend 
on the exact scattering geometry. Figure [2] displays Ra¬ 
man spectra taken under “XX” and “YY” geometries on 


a single-domain crystal. Here, XX denotes that both the 
incoming- and scattered-photon polarizations are along 
the [110]c direction, whereas YY denotes that both po¬ 
larizations are along the [110]c direction; the hexagonal 
c-axis of the crystal is determined to be along the [lll]c 
direction (see below). The measurement configurations 
are illustrated in upper-left inset of Fig. [2j Comparing 
the XX and YY spectra, one sees a clear difference in 
the 380 to 500 cm“^ region, where four peaks can be 
attributed to Ag and Eg modesThe difference is due 
to the fact that both the incoming- and scattered-photon 
polarizations have a non-zero projection along the hexag¬ 
onal c-axis in the XX scattering geometry. This non-zero 
projection leads to a larger (smaller) weight of the Ag 
{Eg) signals in the Raman spectrum compared to that 
obtained in the YY scattering geometry, in which the 
photon polarizations are perpendicular to the hexagonal 
c-axis. 

By performing both XX and YY measurements on only 
one face of the cube and after obtaining the data shown 
in Fig. [21 we can only tell that the hexagonal c-axis is 
along either the [lll]c or the [lll]c direction. In order 
to unambiguously determine the hexagonal c-axis, one 
needs to perform additional measurements on adjacent 
faces of the cube. To verify the validity of our method, 
we have performed six pairs of XX and YY measurements 
on all faces of several single-domain crystals, all yielding 
consistent results supporting the aforementioned picture. 
This gives us confidence on how to identify and orient 
single-domain crystals. Moreover, it suggests the possi¬ 
bility to determine the orientations of individual domains 
in multi-domain samples using Raman spectroscopy. 

For simplicity, here we use the relative inten¬ 
sity ratio between 390 and 425 cm“^, R = 

/(390 cm“^)//(425 cm“^), to represent the key charac¬ 
teristics of the XX and YY Raman spectra in Figs. |2j 
As can be seen from the upper-right inset, R captures 
the most significant difference between the two types of 
spectra: it is approximately 0.8 or 0.5, respectively, when 
the photon polarizations are perpendicular to or partially 
along the hexagonal c-axis. To demonstrate how R can 
be used to characterize surface-domain structures, we 
have performed space-resolved Raman scattering mea¬ 
surements on each of the four types of crystal faces that 
are labeled as A, B, C, and D in Figs. [He-f). The data 
are displayed in Figs.|3l^a2-d2), along with polarized-light 
optical images of the surfaces in the top panels (al-dl). 
In the optical images a horizontal white line indicates the 
trajectory on which the Raman spectra were taken. In 
these Raman measurements, the incident and scattered 
photon polarizations are always kept to be parallel, along 
one of the two perpendicular face diagonals ([110]c and 
[110]c, labeled as “XX” and “YY”, respectively) on the 
surface. 

Indeed, on A-type (Figs. |3l^al-a2)) and D-type 
(Figs. |3l^dl-d2)) faces, R is found to switch between 0.5 
and 0.8 every time the scanning position crosses a bound¬ 
ary between the stripes, and whenever Rxx is around 
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FIG. 2. (Color online) Polarized Raman spectra obtained at 
room temperature, offset for clarity. Data curves are color- 
coded with arrows in the upper-left inset that indicate the 
polarization geometries with respect to the hexagonal c-axis 
(magenta body diagonal of the cube). Vertical dashed lines 
indicate Ag and Eg phonon peaks.— Upper-right inset, XX 
and YY spectra in the 350 to 450 cm“^ region, normalized at 
425 cm“^ 


0.5, Ryy is around 0.8 (and vise versa). The results 
are in perfect agreement with our expectations based on 
the model illustrated in Fig. [TJ On the other hand, R 
is found to remain roughly constant in a given scatter¬ 
ing geometry across the entire B-type (Figs. [3fbl-b2)) 
and C-type (Figs. [3fcl-c2)) faces, for different reasons: 
the B-type surfaces are single-domain, whereas on the C- 
type surfaces the hexagonal c-axes in different domains 
have the same projection onto the surface plane. These 
results render Raman spectroscopy, when used alone, un¬ 
able to distinguish between A-type and D-type, and be¬ 
tween B-type and C-type surfaces. In combination with 
polarized-light microscopy one can easily tell them apart, 
but neither optical inspection nor Raman spectroscopy is 
able to distinguish, e.g.^ between hexagonal c-axis orien¬ 
tations of [lll]c and [lll]c, where the degeneracy is due 
to the fact that light is propagating along the [001] c di¬ 
rection. 

A feasible way to completely determine the domain 
structure by looking at only one surface is to measure 
alternating inclinations, or wrinkles, on the surface. Fig¬ 
ures [3](a3-d3) display surface profiles measured roughly 
along the same trajectory on which we took the Raman 
spectra. We find that the A-type and C-type surfaces ex¬ 
hibit clear zigzag profiles, whereas the B-type and D-type 
surfaces are essentially fiat. Moreover, the angles of the 
zigzag profiles on the A-type (~ 0.80(6)°) and C-type 
1.05(6)°) surfaces are slightly different. To under¬ 
stand these results, we refer to the schematics in Fig. [TJ 


For crystals with {100}c domain walls (Figs.ITj^a, e)), the 
strain tensors of the two domains are Si and S 2 (Eq. [2|); 
when they are contracted with the vector (1, 0, 0), which 
lies within all A-type surfaces and crosses the domains, 
the outcomes of Si and S 2 ((0,(i, d) and (0, —d, —d), 
respectively) have opposite projections along both the 
[010]c and [001]c directions. These are the normal direc¬ 
tions of the A-type surfaces, and thus the surfaces are 
wrinkled. A similar argument can be used to explain 
the wrinkles on C-type surfaces (Figs.[T](b, f)), where the 
strain tensors can be taken as S 2 and S's, which are to be 
contracted with the vector (1, —1,0) and then projected 
along [001]c- For D-type surfaces, we should instead use 
(1,0,0) (or (0,1,0)) as the vector to be contracted with 
S 2 and S's, but here the outcomes have the same pro¬ 
jection along the [010]c (or [100]c) direction, hence the 
D-type surfaces are not wrinkled despite the presence 
of domain walls. No wrinkles are expected on single¬ 
domain B-type surfaces. Using the room-temperature 
lattice constants of CaMn 70 i 2 ^^ the angles of wrinkles 
on A- and C-type surfaces are calculated to be 0.75° and 
1.05°, respectively, which are in good agreement with our 
surface-profile data. 

Taking the above results altogether, we suggest a 
new method to determine ferroelastic domain structures 
based on measurements of only one pseudo-cubic sample 
surface. The most reliable way is to use Raman spec¬ 
troscopy combined with surface profile measurements, 
with polarized-light microscopy being a complementary 
but not necessary method. First, one needs to find the 
orientation of the domain walls by scanning in different 
directions and rotating photon polarizations in the Ra¬ 
man scattering measurement, aiming to maximize the 
contrast in both the Raman and surface profile data. 
Second, the surface domain structure can be know by 
comparing the data with the results shown in Fig. O 
Third, from the surface profile data one can further tell 
apart the aforementioned [lll]c and [llT]c degenerated 
situations on A- and C-type surfaces, using the fact that 
shortened hexagonal c-axis always connects the valleys 
of the wrinkles. Our method is particularly useful when 
the sample is in thin-film form, or when the edges of 
crystals are not along a high-symmetry direction. While 
piezoresponse force microscopy^ is most commonly used 
to study ferroelectric domain structures, our method pro¬ 
vides a route to monitoring the ferroelastic domains both 
above and below the ferroelectric transition temperature, 
which may help improve our understanding of the inter¬ 
play between different ferroic order parameters in multi- 
ferroics. To detect nano-scale domains, tip-enhanced Ra¬ 
man spectroscopy^ii^ can be used to enhance the spatial 
resolution of our method. 


C. Domain switching and detwinning effect 

In addition to methods for characterizing the domains, 
we found that domain structure in CaMn 70 i 2 single crys- 
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FIG. 3. (Color online) (al)-(dl) Optical images of A-, B-, C-, and D-type crystal faces (Figs. [UJe-f)). Horizontal white line 
indicates the trajectory along which the measurements in (a2)-(d2) and (a3)-(d3) are performed. Half-solid-half-dashed line 
indicates the orientation of the hexagonal c-axis (along one of the (lll)c directions) in each domain, with the solid end pointing 
at the top face. (a2)-(d2) Raman intensity ratio R (see text) measured along the trajectories indicated in (al)-(dl), respectively, 
with different photon polarizations. The data are color-coded with the arrows indicating polarization directions in (al)-(dl). 
(a3)-(d3) Surface profiles along the trajectories indicated in (al)-(dl). 


tals can be altered at room temperature under moderate 
uniaxial compression. Figure 0] displays the same face of 
a crystal at different times. Initially, the surface exhibits 
a stripe pattern indicative of presence of {100}c domain 
walls (Fig. m^a)). When a compressive force of about 1.5 
N is applied along the direction shown by the arrows in 
Fig.m^b), the stripe pattern is rotated by 45°, which indi¬ 
cates the formation of new {110}c domain walls and the 
disappearance of the old ones. The magnitude of the ap¬ 
plied force amounts to a uniaxial stress of about 30 MPa 
inside the crystal, and the new stripe pattern persists 
after the force was removed. Similarly, when a compres¬ 
sion is applied along the direction in Fig.lU^c), the stripe 
pattern is found to rotate again by 90°. This domain¬ 
switching behavior can be understood as the following: 
when a uniaxial compression is applied, the ferroelas- 
tic domains rearrange themselves to minimize the length 
along the direction of compression. In the case of com¬ 
pression in the [110]c direction, domains with hexagonal 
c-axes along [lll]c and [llT]c are energetically favored. 
According to the analysis in Section IIII Ai this may re¬ 
sult in the formation of either (OOl)c or (llO)c domain 
walls, but the former would also generate wrinkles on the 
(lOO)c and (OlO)c faces which are incompatible with the 
applied compression. Hence the resulting domain walls 
are parallel to (llO)c. For the same reason, one can rea¬ 


sonably expect that a (lll)c domain can be exclusively 
selected if the compression is primarily along a body di¬ 
agonal of the cube, as is indeed shown to be the case in 
Fig.[l|(d). Apart from taking optical images, we have per¬ 
formed Raman scattering measurements on all six faces 
of the sample in Fig. HJ^d), which consistently show that 
the crystal is highly detwinned. 

This switching behavior of ferroelastic domains in 
CaMn 70 i 2 stems from the fact that the rhombohedral 
distortion is characterized by the shortening rather the 
elongation of a body diagonal, which makes the detwin- 
ning operations simple. Among other rhombohedrally 
distorted perovskites, LaAlOa^ is similar to CaMn 70 i 2 , 
whereas the distortion in BiFeOs features an elongated 
body diagonal, making it rather tricky to prepare certain 
types of domain strncture j^^i^^ 


IV. CONCLUSIONS 

In summary, we have observed and investigated ferroe¬ 
lastic domain structures in single crystals of CaMn 70 i 2 . 
For cube-shaped single crystals with multiple domains, 
we can determine the orientation of individual domains 
by measurements either on two adjacent faces with Ra¬ 
man spectroscopy alone, or on only one face using both 
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FIG. 4. (Color online) Optical images of the same face of 
a crystal: (a) initial state, (b) after compression was applied 
along the face diagonal indicated by the arrows, (c) after com¬ 
pression was applied along the other face diagonal, and (d) 
detwinned state after compressed primarily along a body di¬ 
agonal. 


Raman spectroscopy and stylus surface profiler. The lat¬ 
ter method is also suitable for determining domain struc¬ 
tures in thin-film samples, where only one surface is avail¬ 
able. In addition, polarized-light microscopy provides a 
complementary and convenient way to observe the do¬ 
main structure. Finally, we find that the domain struc¬ 
ture can be altered by moderate uniaxial compression at 
room temperature, which allows for a simple method to 
obtain twin-free samples with a controlled orientation of 
the hexagonal c-axis. Our results offer the opportunity to 
prepare well-defined CaMn 70 i 2 samples, e.g., for spec¬ 
troscopic studies that require single-domain crystals, and 
the methods we use can be readily transferred to stud¬ 
ies of thin-film samples as well as other rhombohedrally 
distorted cubic compounds. 
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